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We study the orbital behavior at the neighborhood of complex unstable periodic orbits in a 3D 
autonomous Hamiltonian system of galactic type. At a transition of a family of periodic orbits 
from stability to complex instability (also known as Hamiltonian Hopf Bifurcation) the four 
eigenvalues of the stable periodic orbits move out of the unit circle. Then the periodic orbits 
become complex unstable. In this paper we first integrate initial conditions close to the ones 
of a complex unstable periodic orbit, which is close to the transition point. Then, we plot the 
consequents of the corresponding orbit in a 4D surface of section. To visualize this surface of 
section we use the method of color and rotation [Patsis and Zachilas 1994]. We find that the 
consequents are contained in 2D "confined tori" . Then, we investigate the structure of the phase 
space in the neighborhood of complex unstable periodic orbits, which are further away from the 
transition point. In these cases we observe clouds of points in the 4D surfaces of section. The 
transition between the two types of orbital behavior is abrupt. 
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1. Introduction 

The aim of this paper is to study the orbital be- 
havior at the neighborhood of a complex unstable 
periodic orbit in a 3D autonomous Hamiltonian sys- 
tem of galactic type. Complex instability is a type of 
instability of periodic orbits that appears in Hamil- 
tonian systems of three or more degrees of freedom. 

In order to study the dynamical behavior at the 
neighborhood of a complex unstable periodic orbit 
we use the method of surfaces of section [Poincare 
1892], which has many applications to Dynamical 
Astronomy [for a review see e.g. Contopoulos 2002]. 
A basic problem in Hamiltonian systems of three de- 
grees of freedom is the visualization of the 4D sur- 
faces of section. Let us assume the phase space of an 
autonomous Hamiltonian system, that has 6 dimen- 
sions, e.g. in Cartesian coordinates, {x,y, z,x,y, z). 
For a given Jacobi constant a trajectory lies on a 
5D manifold. In this manifold the surface of section 
is 4D. This does not allow us to visualize directly 
the surface of section. 

Patsis and Zachilas [1994] proposed a method 
to visualize 4D spaces of section. It is based on rota- 
tion of the 3D projections of the figures in order to 
understand the geometry of the formed structures 
and on color for understanding the distribution of 
the consequents in the 4th dimension. We use for 
this application the "Mathematica" package [Wol- 
fram 1999]. We work in Cartesian coordinates and 
we consider the y = with y > cross section. A set 
of three coordinates (e.g. {x,x,z)) are used for the 
3D projection, while the fourth coordinate (in our 
example z) will determine the color of the conse- 
quents. There is a normalization of the color values 
in the [min(2:), max(2;)] interval, which is mapped 
to [0,1]. Following the intrinsic "Mathematica" sub- 
routines our viewpoint is given in spherical coordi- 
nates. The distance d of the center of the figure from 
the observer is given by "Mathematica" in a special 
scaled coordinate system, in which the longest side 
of the bounding box has length 1. For all figures we 
use d = 1. The method associates the smooth dis- 
tribution or the mixing of colors, with specific types 
of dynamical behavior in the 4th dimension [Patsis 
and Zachilas 1994] [see also Katsanikas and Patsis 
2011]. 

The calculation of the linear stability of a pe- 
riodic orbit follows the method of Broucke [1969] 
and Hadjidemetriou [1975]. We first consider small 
deviations from its initial conditions and then inte- 
grate the orbit again to the next upward intersec- 



tion. In this way a 4D map (Poincare map) is es- 
tablished and relates the initial with the final point. 
The relation of the final deviations of this neigh- 
boring orbit from the periodic one, with the ini- 
tially introduced deviations, can be written in vec- 
tor form as ^ = M ^o- Here ^ is the final devia- 
tion, ^0 is the initial deviation, and M is a 4 x 4 
matrix, called the monodromy matrix. It can be 
shown, that the characteristic equation can be writ- 
ten in the form A"^ + aX^ + PX^ + aA + 1 = 0. Its 
solutions Xi, i = 1,2,3,4, due to the symplectic 
identity of the monodromy matrix, that obey the 
relations A1A3 = 1 and A2A4 = 1 can be written as: 



\y T~ — 



-,z = l,2 



(1) 



where 



hi,2 = 



and 



A = - 4(/? - 2) 



(2) 



(3) 



The quantities 61 and 62 are called the stability 
indices. Following the notation of Contopoulos and 
Magnenat [1985], if A > 0, \hi\ < 2 and I62I < 2, all 
four eigenvalues are complex on the unit circle and 
the periodic orbit is called "stable" (S). If A > 
and > 2, I62I < 2 or < 2, I62I > 2, the pe- 
riodic orbit is called "simple unstable" (U). In this 
case two eigenvalues are on the real axis and two are 
complex on the unit circle. If A > and |6i| > 2 
and 1 62 1 > 2, the periodic orbit is called "double 
unstable" (DU) and the four eigenvalues are on the 
real axis. Finally, if A < the periodic orbit is 
called "complex unstable" (A). In this case the four 
eigenvalues arc complex numbers and they are off 
the unit circle. For the generalization of this kind 
of instability in Hamiltonian systems of N degrees 
of freedom the reader may refer to Skokos [2001]. 
If we have a stable one-parameter (in our system 
the Jacobi constant) family of periodic orbits, the 
four eigenvalues of a stable periodic orbit are com- 
plex on the unit circle. By varying the parameter 
we have a pairwise collision of eigenvalues on two 
conjugate points of the unit circle. From the Krein- 
Moser theorem [e.g. Contopoulos 2002 p. 298] wc can 
decide if after the collision of the eigenvalues they 
will remain on the unit circle and the periodic orbits 
of the family will stay stable, or if the eigenvalues 
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Fig. 1. Stability diagram for —5.5 < Ej < —4.5, that shows 
the stabiUty of family xl (with black line) and its bifurcating 
family of p.o. xlvl (with red line). 



will move out from the unit circle into the com- 
plex plane forming a complex quadruplet. In this 
latter case the periodic orbits of the family will be- 
come complex unstable and we will have a transition 
from stability to complex instability (also known as 
Hamiltonian Hopf Bifurcation) . From an analytical 
point of view, the transition to complex instabil- 
ity has been studied using the Hamiltonian itself 
[Heggie 1985, Broer et al. 2007, Olle et al. 2008] or 
4D symplectic maps [Bridges et al. 1995]. In both 
cases, the approach consists of normal forms tech- 
niques [Olle et al. 2005a, b] to simplify the Hamil- 
tonian (or the map) and describe the local phase 
space structure near the critical periodic orbit (or 
fixed point in the discrete context). Such analysis 
shows, that the transition to complex instability 
gives rise to bifurcating invariant 2D tori in the flow 
context or invariant curves in 4D symplectic maps 
(as Poincare map). The numerical computation of 
these invariant objects has been done for Hamil- 
tonian systems of three degrees of freedom [Pfen- 
niger 1985b, Olle and Pfenniger 1998, Olle et al. 
2004] and for 4D symplectic maps [Pfenniger 1985a, 
Jorba and Olle 2004]. It is remarkable, that there 
exists not one but two kinds of Hamiltonian Hopf 
bifurcations (as it happens in the usual dissipative 
setting), depending on the coefficients of the nor- 
mal form [Van der Meer 1985]. These two kinds of 
bifurcations are usually called direct (supercritical) 
and inverse (subcritical). From a numerical point of 
view we can distinguish the two kinds of Hamilto- 
nian Hopf bifurcation if we take firstly initial condi- 
tions in the vicinity of a complex unstable periodic 
orbit near the transition point from stability to com- 
plex instability. After that, we plot the consequents 
of the corresponding orbit in the surface of section. 



If the consequents are confined, we have a direct 
Hamiltonian Hopf Bifurcation. Otherwise, if the or- 
bit escapes we have an inverse Hamiltonian Hopf 
Bifurcation. In this paper we examine the structure 
of the invariant surfaces in the 4D surface of section 
in the neighborhood of a complex unstable periodic 
orbit after a direct Hamiltonian Hopf bifurcation. 
The inverse Hopf bifurcation is not considered in 
the present paper, because in such a case all orbits 
close to a complex unstable periodic orbit escape 
[Jorba and Olle 2004] and no confined structures 
appear. 

2. The Hamiltonian System 

The system we use for our applications is described 
in details in Katsanikas and Patsis [2011]. It ro- 
tates around its z-axis with angular velocity rif, = 
60 km kpc~^. The Hamiltonian of our system 
is: 

H{x,y,z,x,y,z) = 
^{x^ + f + i'^) + ^{x,y,z) 

-lnl{x' + y% (4) 

where $(x,y,z) is our potential. The potential 
^{x,y,z) in its axisymmetric form can be consid- 
ered as a representation of the Milky Way approxi- 
mated by two Miyamoto disks with masses Mi and 
M2 respectively [Miyamoto and Nagai 1975]. The 
parameters we use in the present study are the same 
as in Katsanikas and Patsis [2011]. In our units, 
the distance R=l corresponds to 1 kpc. The veloc- 
ity unit corresponds to 209.64 km/ sec. For the Ja- 
cobi constant (hereafter called the "energy") Ej=l 
corresponds to 43950 (km/sec)'^, while the velocity 
unit is 209.64 km/ sec. 

3. Spaces of section 

A method to follow the stability of a family of peri- 
odic orbits in a system is by means of the "stability 
diagram" [Contopoulos and Barbanis 1985, Pfen- 
niger 1985a]. The stability diagram gives the evolu- 
tion of the stability of a family of periodic orbits in 
a system as one parameter varies, by means of the 
evolution of the stability indices 61 , 62 . In our case 
the parameter that varies is the energy Ej. Fig. [T] 
gives the evolution of the stability of the central 
family of periodic orbits xl on the equatorial plane 
[Contopoulos and Papayannopoulos 1980], and its 
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Fig. 2. 2D projection (x,z) of the 4D surface of section at 
the neighborhood of a complex unstable (A) periodic orbit 
of the 3D family xlvl for Ej — —4.619 (for 70 intersections). 



bifurcations for —5.5 < Ej < —4.5. We observe, 
that xl is initially stable and at Ej= —5.1644 it 
becomes simple unstable. There we have a 5 — >• C/ 
transition and a new 3D family, xlvl [Skokos et al 
2002a,b], is bifurcating and is stable. For £^j=— 4.62 
the family xlvl becomes complex unstable and the 
stability indices meet each other. After that the val- 
ues of bi^2 are complex numbers. 

We investigate the structure of the phase space 
by perturbing the initial conditions of a complex 
unstable periodic orbit in the 4D surface of section 
after a transition from stability to complex instabil- 
ity. For example we can apply a perturbation Ax = 
10~^ at the initial conditions of a complex unsta- 
ble periodic orbit for Ej=— 4.619 with initial condi- 
tions (xo,xo,^o,^o) = (0.275137727,0,0.359838,0). 
Then, the x initial condition will be xq + Ax. From 
previous papers [Contopoulos et al. 1994, Papadaki 
et al. 1995] we know that, for a number of inter- 
sections, a spiral appears in the consequents at the 
neighborhood of a complex unstable (A) periodic 
orbit in 2D projections of the 4D surfaces of sec- 
tion. This structure is encountered also in the case 
we study (Fig. [2|). We observe that the consequents 
are confined and this means that we have a direct 
Hamiltonian Hopf bifurcation. 

We apply the method of color and rotation to 
study the invariant surface in the neighborhood of 
the complex unstable periodic orbit. In Fig. [3] we 
present the orbital behavior close to the complex 
unstable periodic orbit for 50 intersections with the 



surface of section. We use the 3D 



^x, X, z ) projec- 



tion for plotting the points and the i values to 
color them. Our point of view is (9, (p) = (32°, 53°). 
We observe a spiral (with three spiral arms) in the 
neighborhood of the complex unstable periodic or- 




Fig. 3. The orbital behavior close to a complex unstable 
(A) periodic orbit of the 3D family xlvl for Ax = 10"'* and 
Ej = —4.619 (for 50 intersections). We use the {x, x, z) space 
for plotting the points and the z value to color them. Our 
point of view in spherical coordinates is given by {9, (p) = 
(32°, 53°). 




0.06 0.31 

Fig. 4. The same case as Fig. |3] for 2000 intersections. 



bit. The first 13 points are very close to the center 
of the spiral and are not discernible in the figure. 
The 14th point is located on a spiral arm, the 15th 
point is located on another spiral arm and the 16th 
point is located on the third spiral arm. The 17th 
point is on the first spiral arm again and so on. In 
this way the points fill the space. Along every spiral 
arm we have a succession of colors (Fig. [3|). In Fig. [3] 
we observe that if we begin from the center along 
the spiral arm going through the 15th point we have 
a succession of colors from green to light blue, light 
blue to blue and blue to violet. Along the spiral arm 
through the 16th point we observe a smooth color 
variation from green to light blue, to blue, to light 
blue, to green, to green-yellow, to orange-red. Fi- 
nally along the third spiral arm we observe also a 
smooth color variation. Starting from the center of 
this spiral arm, we see that green becomes green- 
yellow, then yellow, green-yellow, green and finally 
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light blue (Fig. [3]). 

Figure H] depicts the same case as in Fig. [3] 
and we observe the orbital behavior close to the 
complex unstable periodic orbit for 2000 intersec- 
tions with the surface of section. We cannot distin- 
guish the spiral any more, but we observe a disk- 
like structure. The central area is saturated due to 
the presence of many points. This object has a very 
small thickness and is also slightly warped. Here- 
after we will refer to this structure as the "disk" 
or "disk structure". In the terminology of Jorba & 
Olle [2004] such structures are called confined tori. 
Small deviations from the pure planar geometry are 
due to a warp of the disk in the 3D space. On this 
disk structure, from left to right, we have a suc- 
cession of colors from red to orange, orange to yel- 
low, yellow to green, green to light blue, light blue 
to blue and blue to violet. For 10000 intersections 
we observe that the points fill the disk structure 
(Fig. [5]), which has a smooth color variation. The 
smallest distance of the consequents on this disk 
structure from the complex unstable periodic orbit 
is 0.000042 and the largest distance is 0.059764 in 
the 3D projection (x, x, z) of the 4D surface of sec- 
tion. The consequents reach an outermost distance, 
then they move inwards reaching a minimum dis- 
tance, then they move outwards etc. Some values 
for the minimum and maximum distances are given 
in Tables 1 and 2. We underline the fact, that the 
disk has an internal three-armed spiral structure. 
This means that if we choose arbitrarily a point on 
the disk, the subsequent consequents will follow a 
spiral pattern as the one presented in Figs. El |4] and 
m which lies on the disk. For this orbit we calcu- 
lated also the "finite time" Lyapunov Characteristic 
Number (LCN), i.e. 



LCN{t) = J In 



m 



e(to 



where ■^(io) and (,{t) are the distances between two 
points of two nearby orbits at times t = and t 
respectively( see e.g. Skokos [2010]). We found, that 
initially it decreases and finally it levels off around 
6 X 10~^ after about 4500 intersections (Fig. [6l lower 
blue curve). 

We observe the same behavior (the same disk 
structure) if we apply a perturbation Ax = 10~^ 
to the initial conditions of the complex unstable 
periodic orbit for Ej=— 4.619 in the same way we 
did in the previous case for Ax = 10~^. The only 
difference we find is that the minimum distance of 
the consequents from the complex unstable periodic 



0.24 



0.31 




0.375 



345 



0.08 



Fig. 5. The same case as Fig. [3] for 10000 intersections. 




1x10* 



t 



Fig. 6. The evolution of the two LCN{t). The blue curve 
corresponds to the case of the confined torus, while the red 
to the case of the cloud. The axes are in logarithmic scale. 



0.38, 




Fig. 7. The orbital behavior close to a complex unstable 
(A) periodic orbit of the 3D family xlvl for Ax = 10^^ and 



-4.619 (for 10000 intersections) 
spherical coordinates is given by (6, cj>] 



Our point of view in 
= (15°, 90°). 
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Table 1. Some examples of points on the disk for a perturba- 
tion Aa; — 10"** from the complex unstable periodic orbit for 
Ej — —4.619, which have maximum distance from the periodic orbit. 



nth point 


221 


522 


802 


2679 


8904 


max distance 


0.056284 


0.058844 


0.058453 


0.059693 


0.059764 



Table 2. Some examples of points on the disk for a perturba- 
tion Ax = 10^** from the complex unstable periodic orbit for 
Ej = —4.619, which have minimum distance from the periodic orbit. 



nth point 


300 


600 


912 


2780 


8990 


min distance 


0.000338 


0.000316 


0.000042 


0.000050 


0.000173 




0.5 



Fig. 8. The orbital behavior close to a complex unstable 
(A) periodic orbit of the 3D family xlvl for Ax = 10~^ and 
Ej = —4.619 (for 310 intersections). Our point of view in 
spherical coordinates is given by (6^,0) = (15°, 90°). 




X 



Fig. 9. The 3D projection {x,x,z) of the consequents close 
to a complex unstable (A) periodic orbit of the 3D family 
xlvl for Aa; = 10~^ and Ej = -4.619 (for 3000 intersec- 
tions). Our point of view in spherical coordinates is given by 
(e,<^) = (55°,95°). 



orbit is now 0.000999 and the maximum distance 
is 0.0592138. The minimum distance is 24 times 
larger than the minimum distance for a perturba- 
tion Ax = 10"*^ and for this reason we have a small 
hole at the center of the disk structure. Now if the 
perturbation is even larger, Ax = 10~^, we observe 
a different disk structure (Fig. [7]). The new disk be- 
comes thicker than the previous disk we have stud- 
ied with initial Ax = 10^^ or Ax = 10~^ away from 
the periodic orbit. Now the "thickness" in the three 
coordinates is (AX,AX,AZ) = (0.09,0.2,0.04) 
instead of {AX,AX,AZ) = (0.07,0.16,0.03) we 
had before. On this "disk structure" we observe 
a smooth color variation from red to violet (Fig. 
[7|. This means that the 4th dimension follows the 
topology of this structure in the 4D surface of sec- 
tion. The smallest distance of the consequents on 
this structure from the complex unstable periodic 
orbit is 0.009999 and the largest distance is 0.078655 



in the 3D projection (x, x, z) of the 4D surface of 
section. Thus, the smallest distance of the conse- 
quents on the disk structure for the perturbation 
Ax = 10^^ is now 238 times larger than the small- 
est distance of the Ax = 10~^ perturbation. As ex- 
pected the hole in middle of the "disk" is now even 
larger (Fig. [7]) . 

Now we increase the perturbation by tak- 
ing Ax = 10-\ always for =-4.619. The 
first approximately N=300 consequents form a 
toroidal surface. During this period the "finite time" 
LCN{t) of the orbit decreases to a value 1.75 x 10~^ 
(Fig. m upper, red curve). Beyond that point it fluc- 
tuates as the time increases and finally increases 
and tends to level off around 2.5 x 10^^ after about 
4700 intersections. On our 4D surface of section the 
consequents fill a cloud around the toroidal object 
during the same time. 
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Fig. 10. Diagram of the AZ variation of the consequents 
close to a complex unstable periodic orbit versus Ej. 



In Fig. [8] we observe the first 310 consequents. 
Besides the points that form the toroidal object, we 
have included also the very first points that depart 
from it. These points correspond to the cuspy fea- 
tures we can observe in the red dense ring in Fig. [9l 
In Fig. [8] the 4th dimension z is represented by the 
colors. Smooth color variation in the consequents is 
again present and this means that we have a toroidal 
surface in the 4D space of section. By considering 
more points on the surface of section (Fig. [9]), we 
observe that they deviate from this toroidal surface 
and they soon occupy a large volume of the phase 
space. This indicates stickiness [Contopoulos and 
Harsoula 2008]. 

As mentioned previously, the maximal LCN for 
the case of the cloud (Fig. ^ is around 2.5 x 10~^, 
while the same index for the orbit of the confined 
torus (Fig. [5]) is lower, 6 x 10~^. This difference is 
expected by the nature of the two orbits. In the 
cloud we have scattered points in 4 dimensions, 
while in the case of the confined torus, chaos is due 
to the different maximum and minimum distances 
from the center reached as the orbit "fills" the disky 
structure of this torus (see Tables 1 and 2). 

We choose now a Ax = 10~^ deviation from the 
initial conditions of the periodic orbit, that gives 




X 



Fig. 11. The 3D projection (a;, i, z) of the consequents close 
to a complex unstable (A) periodic orbit of the 3D family 
xlvl for Ax = 10"'' and Ej = -4.50 (for 4000 intersections). 
The periodic orbit is depicted by a black point. Our point of 
view in spherical coordinates is given by (S, 0) = (62°, 5°). 




Fig. 12. The orbital behavior close to a complex unstable 
(A) periodic orbit of the 3D family xlvl for As = 10~^ and 



Ej = -4.50 (for 2000 intersections), 
spherical coordinates is given by (6, cj> 



Our point of view in 
= (55°, 88°). 



a confined torus at 



-4.619 and we increase 



Ej. In Fig. [TD] we depict the variation of the thick- ues of the energy larger than —4.59. For example for 
ness AZ of the consequents in the 4D surface of sec- 
tion versus the energy Ej. We observe that for val- 
4.62 and —4.59 we have small 



ues of Ej between 
values of AZ and for values of Ej between —4.59 
and —4.42 we have large values of AZ. The first 
interval of Ej corresponds to the disk structures 
that we described before. The second interval corre- 
sponds to the clouds of points that we find for val- 



Ej = —4.50 we observe a cloud of points with values 
of z from —1.5 to 1.5 (Fig. Ill|) . The periodic orbit 
is depicted black point and it has initial condi- 
tions (xo,i:o,-zo,io) = (0.307260,0,0.432148,0). In 
Fig. [12] we observe a mixing of colors for this cloud 
and this indicates a chaotic behavior in 4D space. 

Let us investigate now the changes observed in 
the neighborhood of the complex unstable periodic 
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Fig. 13. The orbital behavior close to a complex unstable 
(A) periodic orbit of the 3D family xlvl for Ax = 10"'* and 
Ej = —4.59 (for 1250 intersections). Our point of view in 
spherical coordinates is given by {0,(j)) = (15°, 90°). 



1.5 



0.5 




Fig. 14. The 3D projection (x, x, z) of the consequents close 
to a complex unstable (A) periodic orbit of the 3D family 
xlvl for Ax = 10"'' and Ej = -4.59 (for 4000 intersections). 
Our point of view in spherical coordinates is given by (6, (p) = 
(55°, 88°). 



orbit when we vary the energy. At the transition 
point of Fig. [TOl for Ej = —4.59 and for a pertur- 
bation of the initial conditions equal to Ax = 10"^, 
the first 1250 consequents remain on a disk struc- 
ture as we can see in Fig. [131 We note that the points 
do not suffice to fill densely the "disk" . On this disk 
structure we observe a color succession from violet 
to blue, to light blue, to green, to yellow, to or- 
ange, to red. We underline that the succession of 
colors means that the 4th dimension z of the con- 
sequents lies on this disk structure. In Fig. [T3] we 
see that if we integrate the orbit for more time, the 
consequents start to leave this disk structure (that 
is depicted with black color) and form a cloud of 
points around it. This is again a typical behavior of 
a sticky orbit. 

4. Conclusions 

In this paper we studied the phase space structure 
at the neighborhood of a complex unstable periodic 
orbit after a direct Hamiltonian Hopf bifurcation. 
Close to a Hamiltonian Hopf bifurcation we observe 
that: 

(1) We have a spiral structure formed by the conse- 
quents at the neighborhood of complex unstable 
periodic orbits in the 2D projections of the 4D 
surface of section as in Contopoulos et al [1994] 
and Papadaki et al [1995]. Here we find that 
we have this spiral structure in the 3D projec- 
tions and we observe a smooth color variation 



along their arms. This means that this spiral 
structure is a 4D object. 

(2) The consequents near a complex unstable peri- 
odic orbit arrive at a maximum distance from 
the periodic orbit and they move inwards. 
By repeating this process, a disk structure is 
formed in the 3D projection of the 4D surface of 
section, which is called a confined torus [Pfen- 
niger 1985a,b, Jorba & Olle 2004 and Olle et 
al 2004]. On this disk structure we observe a 
smooth color succession. This means that we 
have a disk structure (the confined torus) also 
in the 4D space of section. 

(3) If we apply larger perturbations to the initial 
conditions we observe, that the disk structures 
become toroidals with smooth color variation 
and holes at their centers. For a critical value of 
the Ax perturbation of the initial conditions we 
observe for a number of intersections a toroidal 
surface with smooth color variation. However, 
later, the consequents leave this toroidal sur- 
face and move away occupying a larger volume 
in the phase space. This is a case of stickiness 
[Contopoulos & Harsoula 2008] in the case of 
complex instability and confined tori. This is 
the first time that we visualize stickiness in the 
neighborhood of a complex unstable periodic 
orbit. 

(4) As the value of energy increases we do not find 
confined tori anymore, but clouds of points in 
the neighborhood of complex unstable periodic 
orbits. These clouds have mixing of colors. This 
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means that we have strong chaos in the 4D 
space of section when we go far from the transi- 
tion point from stabihty to complex instabihty. 
(5) The calculation of the "finite time" LCN, gives 
a global value of the variation of the volume 
filled by the consequents of the orbits. We found 
that the LCN{t) curve levels off after more 
than 4500 intersections and is around 2.5 x 10^^ 
for the case of the cloud (Fig. [9|, and around 
6 X 10~^ for the confined torus (Fig. [5]). On the 
other hand, the method of color and rotation 
gives in a much more direct and detailed way 
the volumes of the phase space occupied by an 
orbit during its integration. This is of particular 
importance in Galactic Dynamics. 
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